Recollisions and correlated double ionization with circularly polarized light 
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It is generally believed that the recollision mechanism of atomic nonsequential double ionization is 
suppressed in circularly polarized laser fields because the returning electron is unlikely to encounter 
the core. On the contrary, we find that recollision can and does significantly enhance double ion- 
ization, even to the extent of forming a "knee", the signature of the nonsequential process. Using 
a classical model, we explain two apparently contradictory experiments, the absence of a knee for 
helium and its presence for magnesium. 

PACS numbers: 32.80.Rm, 05. 45. Ac 



Multiple ionization of atoms and molecules is usually 
treated as a rapid sequence of isolated events. However, 
in the early 90's, experiments using intense laser pulses 
found double ionization yields which departed from these 
treatments by several orders of magnitude [l|, thereby 
casting doubt on the uniqueness of the sequential (uncor- 
rclatcd) multiple ionization channel. These results con- 
stitute one of the most striking surprises of recent years in 
intense laser-matter interactions, and this new ionization 
channel ("nonsequential double ionization" (NSDI) [l[) 
has emerged as one of the most dramatic manifestations 
of electron-electron correlation in nature ■ 

Most of the experimental observations of this striking 
process [l[ used linearly polarized (LP) laser fields. In 
this setting, the precise mechanism that makes electron- 
electron correlation so effective follows the recollision (or 
"three-step" ) scenario 0, H| : An ionized electron, after 
picking up energy from the field, is hurled back at the ion 
core upon reversal of the field and dislodges the second 
electron. NSDI has become an integral part of attosec- 
ond physics [j| since the recollision mechanism requires 
high-intensity, short-pulse lasers. In fact, it is hoped that 
attosecond control will provide insights into the dynam- 
ics of electron-electron collisions as well as complex mul- 
tielectron collision phenomena. 

Many questions remain unanswered regarding strong- 
field double ionization, and one that is still completely 
open concerns polarization. The stakes are high when 
it comes to understanding the influence of polarization 
since it is well known that the emission of harmonics in 
atoms and molecules is strongly dependent on the ellip- 
ticity of the driving field [H, 0] , which can therefore act as 
a control knob. Indeed, recent work shows that elliptic 
polarization (EP) provides a new control mechanism for 
recollision physics in high- field ionization Q. To illus- 
trate the conceptual difficulty associated with polariza- 
tion, consider circularly polarized (CP) fields: The rec- 
ollision scenario which works so well in LP fields is much 
more difficult to justify in fields in which the ionized elec- 
trons tend to spiral out from the core and to miss it 
Therefore one has to come to expect any recollision in CP 
fields to be due to contamination by small amounts of EP 
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FIG. 1: Computed double ionization probabilites for Hamil- 
tonian {1} as a function of the intensity of the field I for 
ui — 0.0584 a.u. (780 nm) and circular polarization, a = 3 
and b — 1 (magnesium case). The inset displays the same 
computation for helium (a = b — 1). 



fields, but not otherwise. The matter would rest there 
if it were not for conflicting experimental evidence: In 
some experiments using CP fields, the double ionization 
yields follow the sequential mechanism, confirming cur- 
rent thinking, whereas in others these yields are clearly 
several orders of magnitude higher than expected, in ap- 
parent contradiction with it. The question we resolve 
here is: Are recollisions possible in pure CP fields or does 
one have to rely on a small residual ellipticity? 

In order to benchmark the process, multiple ioniza- 
tion yields are computed as a function of the intensity 
of the laser. The signature of the nonsequential double 
ionization is a characteristic "knee" shape in the proba- 
bility versus intensity curve (see Fig.Q]). Experimentally, 
a knee has been reported in for nitrogen oxide [l(| and 
for magnesium pdl |. Other experiments (including those 
on helium fi~2| ]) conclude that no knee is formed in CP 
fields. 

In what follows we explain these seemingly contradic- 
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tory findings using a classical Hamiltonian model. In 
addition, we show that, contrary to common belief, rec- 
ollision can be the dominant mechanism leading to en- 
hanced double ionization yields. The presence or absence 
of a knee amounts to finding the conditions for the rec- 
ollision mechanism to work in CP fields. It turns out 
that enhanced double ionization in CP fields belongs to 
a class of rotational problems which are closely related to 
ionization of Rydberg states in CP microwave fields 13 1 
as well as asteroid capture [l4| and matter transport 15 j 
in celestial mechanics. 

The CP field is special among EP fields because viewed 
from a frame rotating with it, it turns into a con- 
stant static field, and the problem becomes analogous 
to the combined Stark-Zeeman effect (or the so-called 
crossed fields problem [lj| [lj}) with a substantial Cori- 
olis term which leads to subtleties where the energy is 
concerned [IH, . The key aspect in this frame change 
is the emergence of the Stark saddle [III [2(| over which 
pre-ionized electrons penetrate the core or excited core 
electrons move away from it. Under special conditions 
the ionized electron can reach the core and collide with 
the inner electron, leading to an appreciable amount of 
NSDI and its characteristic knee signature. 

Entirely classical interactions turn out to be adequate 
to generate the strong two-electron correlation needed 
for double ionization [2|, |21| . Therefore we work with 
the classical Hamiltonian model of pseudo-two electron 
atoms with soft Coulomb potentials [22j]. The Hamilto- 
nian is given by (23|: 



n = 



|Pi| 2 + |P2| 2 



xi - x 2 
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xi| + y|x 2 | 
+ (xi+x 2 )-E(i), 



(1) 



where x^ = (xi,yi) is the position of the i-th elec- 
tron, p s ; = {p x ,i,Py,i) is its canonically conjugate mo- 
mentum and |-| denotes the Euclidean norm in K 2 . 
The laser field is circularly polarized, i.e. E(i) = 
Eof (t) (e x sin cot + e y cos cot) where Eq is the ampli- 
tude of the laser field, cu its frequency chosen as to = 
0.0584 a.u.(780 nm), and / (t) the envelope of the pulse 
with a two laser cycle ramp-up, six laser cycle plateau 
and two laser cycle ramp down. The two parameters a 
and b used in the soft Coulomb potential are determined 
such that b is sufficiently small to allow significant energy 
exchange during recollisions, and a such that in the ab- 
sence of the field, all the initial conditions are bounded 
(to prevent self- ionization) . Therefore the choice of a is 
related to the binding energy of the ground state: For 
He suitable parameters are a = 1 and 6=1 for a ground 
state energy E g = —2.24 a.u. (24|. For Mg, the ground 
state energy is £ g = —0.83 a.u., and we choose a = 3 and 
b = 1 . We consider a large assembly of initial conditions 
(typically 1.2 x 10 5 ) in the microcanonical ensemble of 
Hamiltonian ([TJ) in the absence of the field [25[. A sta- 




FIG. 2: NSDI in magnesium with a laser intensity I = 2 x 
10 13 W • cm . The upper panel displays the energy (kinetic 
energy plus Coulomb interaction with the nucleus) of each 
electron (red and blue curves) as a function of time and the 
green curve represents the energy of the (Coulomb) repulsion 
between the two electrons. The arrow indicates the moment 
of recollision. Lower left panel: The position of each electron 
in the rotating frame. Lower right panel: The position of 
each electron in the rotating frame during the time interval 
indicated by the dashed lines in the upper panel. The arrows 
indicate the positions of the electrons at the recollision and 
the direction of their motion. In the lower panels we display 
the position of the saddle point by a cross. 



tistical analysis of the trajectories provides the double 
ionization yield versus intensity. The results, reported in 
Fig. [I] show a knee for Mg, the characteristic features 
of which are in very good agreement with experimen- 
tal findings [llj : The correlated double ionization yield 
peaks at approximately 3 x 10 13 W • cm~ 2 . A similar cal- 
culation performed for He (inset of Fig. [1]) does not give 
any evidence of a knee, once again in agreement with 
experiment (l2| . 

For Mg, with a CP laser field intensity I = 2x 10 13 W- 
cm~ 2 it turns out that approximately 6% of the trajecto- 
ries are subjected to a strong recollision (i.e., an electron 
is ionized and then returns to approach the other one 
closer than 3.2 a.u.). Thus, in this case, counterintu- 
itively, recollisions are not marginal and do play a lead- 
ing role in the observed knee. In what follows, we use 
tools from nonlinear dynamics [25[ to identify the mech- 
anism by which the electrons undergo correlated double 
ionization. In a nutshell, the correlated double ionizing 
trajectories follow a modified three-step model: First, an 
electron is ionized by escaping through a saddle. If it is 
trapped close to a periodic orbit (this family of outer- 
electron periodic orbits will be explained in the next sec- 
tion), it might return to the core through the saddle, 
rccollide with the core electron, and result in correlated 
double ionization when both electrons escape, once again, 
through the saddle. In Fig. [5] we display a typical exam- 
ple of such a trajectory. 

We have noted before [25[ that, in the absence of the 
field, a typical two-electron trajectory is composed of one 
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electron close to the nucleus (the "inner" electron) and 
another further away (the "outer" electron), with quick 
exchanges of the roles of each electron. This distinction 
is crucial when the laser field is turned on : The outer 
electron is mainly driven by the field, whereas the field 
competes with the nuclear interaction in driving the inner 
electron. Therefore we can use reduced models to identify 
the route the outer electron follows to the core. As in 
the LP case [25[ , the cooperation of the inner and outer 
electrons is essential for ionization in CP. 

Outer electron dynamics- As for LP, finding recollid- 
ing trajectories amounts to finding the initial conditions 
xo and po in phase space such that the trajectory re- 
turns to its initial position (in the core region) after 
some recollision time t TOCO n, i.e., by solving the system 
x(irecoii) = xo. A reduced model for the outer electron 
dynamics is obtained by neglecting the interaction with 
the nucleus and with the other electron (25[. In this case, 
the maximum recolliding kinetic energy brought back by 
the outer electron is nU p where k » 3.17 [1,111 (and U p 
is the ponderomotive energy) as it is exactly the case for 
LP. The analysis of recolliding trajectories shows that 
this is not the dominant behavior observed for the outer 
electron and therefore is not the main mechanism for rec- 
ollisions. The correct mechanism emerges from viewing 
the reduced model in a rotating frame, where we compute 
the distance p(t) between the position of the electron and 
(E /u 2 , 0) as 

p(t) = ||x + p t\\, 

where x = x — (E /uj 2 , 0) and p = Po — (0, Eq/u). For 
Po = 0, the outer electron follows a circular periodic orbit 
(with constant momentum) and so stays in the vicinity 
of the core region. This is actually what is seen from a 
typical recolliding trajectory: po ~ so that it is located 
in the neighborhood of one of the periodic orbits (with 
p(t) constant, see lower left panel of Fig. [2J. The con- 
dition po ~ gives a natural criterion for the feasibility 
of recollisions: If the condition corresponds to admissible 
initial conditions, recollisions are possible. In contrast, if 
the condition po ~ is far away from admissible initial 
conditions, no recolliding trajectories are possible. For 
780 nm and I — 2 x 10 13 W • cm~ 2 the recollision condi- 
tion gives po ~ 0.4 while for / = 10 15 , it gives po ~ 2.9. 
The first case corresponds to admissible conditions and 
thus explains why we observe recollisions in Mg. In con- 
trast, the second case is not allowed which explains the 
absence of recolliding trajectories for He. 

Starting close to the circular periodic orbit, the outer 
electron finds its way through the saddle point to the 
core region because of the Coulomb interaction with the 
nucleus, and might collide with the other electron. De- 
pending on the efficiency of the collision [26[, both elec- 
trons may leave the core, resulting in NSDI (see lower 
right panel of Fig. [2]) . A model for the outer electron dy- 
namics which takes into account the interaction with the 
nucleus is obtained from Hamiltonian (JT|) by neglecting 
the interaction with the other electron. This dynamics is 




FIG. 3: Poincare sections of the inner electron reduced model 
in the 3D space (x, y, P^> = xp y — yp x )- The admissible set is 
represented by the surface and Poincare sections are repre- 
sented by dots. Lower panel: Sections for a Jacobi constant 
K, = —0.659 which corresponds to the constant for the inner 
electron in Fig. [2]when it is close to the nucleus and the outer 
electron is far away. For comparison, we display the Poincare 
section associated with this part of the trajectory for the full 
model (crosses). Upper panel: Poincare surface for a Jacobi 
constant K, = —0.5, when the ionization channel is open. The 
arrow represents the jump, from an invariant torus up to the 
chaotic unbounded area, experienced by the inner electron 
during NSDI. 

best appreciated in the rotating frame where the reduced 
Hamiltonian becomes 

Ip! 2 

1C=— - uj(xp y - yp x ) + E x, (2) 

V>l 2 + « 2 

where x = (x,y) is the position of the electron in the 
rotating frame, p = {p x ,p y ) its canonically conjugate 
momentum, and K is the Jacobi constant [27| . The 
saddle point is located at i = a:*, y — 0, p x = 
and p y = u>x* where is a real (negative) solution of 
lo 2 x — E — 2x/(x 2 + a 2 ) 3 / 2 = which corresponds to a 
saddle of the so-called zero- velocity surface [27|. When 
the ionization channel is open (depending on the Jacobi 
constant, see Fig. [3]), it becomes possible for an elec- 
tron to penetrate the core region. This is exactly what 
happens to the outer electron which is in the vicinity of 
the circular periodic orbits. The trajectory displayed in 
Fig. [2] (lower left panel) shows that the outer electron 
enters the core region close to the saddle indicated by a 
cross. 

Inner electron dynamics- Hamiltonian @ is also a 
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model for the dynamics of the inner electron, but it is 
only valid in the core region where the effect of the po- 
tential is strong so that the inner electron is bound. In 
order to investigate its dynamics, we construct Poincare 
sections of Hamiltonian ([2]) with a surface of equation 
xp x + UPy = (see Fig. [3j. For the relevant range of Ja- 
cobi constants, the dynamics shows an elliptic island in 
the core region which binds the inner electron. At some 
critical Jacobi constant an ionization channel opens up, 
raising the possibility for the inner electron to leave this 
bound region after a small exchange of energy during a 
recollision with the outer electron. The ionization chan- 
nel through a saddle is similar to the one observed for 
the ionization of Rydberg atoms in strong CP microwave 
fields through electron collisions with the core [l3j |. The 
behavior of the inner electron is determined by its Ja- 
cobi constant: For small Jacobi constants, the accessible 
set is composed of two disconnected surfaces, a bounded 
one close to the core and another, unbounded one which 
corresponds to ionized positions. Because of the choice 
of initial conditions for the atom - on the ground state 
energy surface - this second component is initially not ac- 
cessible to the inner electron. However, for larger Jacobi 
constants, the two surfaces merge into a single one, thus 
establishing a connection between the core and unbound 
regions. 

When the Jacobi constant is small, the set which is 
accessible to the inner electron is organized by invariant 
tori which fill up all the accessible domain of the inner 
electron (see Fig. [3J lower panel). But when the con- 
stant increases, some invariant tori break up, leading to 



coexistence between regular structures and a chaotic sea 
(see upper panel, Fig. [3]). This chaotic sea, connected 
to the unbound region, defines an ionization channel for 
the inner electron. This channel passes close to a saddle 
point which gives a natural guide for the ionization seen 
in the rotating frame. The way for the inner electron to 
move from an invariant torus associated with a small Ja- 
cobi constant up to the chaotic sea is through collisions 
with the outer electron when it returns to the core re- 
gion. This jump in the Jacobi constant is symbolized by 
an arrow in Fig. [3J The higher the Jacobi constant is, 
the larger is the channel of ionization around the saddle 
point. Thus, when the channel is very narrow, the inner 
electron may be trapped for a while before finding its way 
out and ionize, leading to a delay between the last colli- 
sion and ionization of the inner electron. The amount of 
energy the inner electron gains from the returning elec- 
tron determines its chances to ionize by leaving an in- 
variant torus to reach the chaotic sea where it is swept 
away by the laser field (see Fig. [3]). Recollisions lead 
to double ionization if the inner electron gains enough 
energy to jump into the chaotic (unbound) region while 
leaving the outer electron with enough energy to remain 
ionized. Such special recollisions are the main source of 
enhancement of the double ionization probability with a 
circularly polarized field and the formation of a knee with 
varying intensity. 
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